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Abstract We performed for the first time stereoscopic triangulation of coronal 
loops in active regions over the entire range of spacecraft separation angles 
{asep ^ 6°, 43°, 89°, 127°, and 170°). The accuracy of stereoscopic correlation 
depends mostly on the viewing angle with respect to the solar surface for each 
spacecraft, which affects the stereoscopic correspondence identification of loops 
in image pairs. From a simple theoretical model we predict an optimum range 
of asep ^ 22° — 125°, which is also experimentally confirmed. The best accu- 
racy is generally obtained when an active region passes the central meridian 
(viewed from Earth), which yields a symmetric view for both STEREO space- 
craft and causes minimum horizontal foreshortening. For the extended angular 
range of agep ~ 6° — 127° we find a mean 3D misalignment angle of jj^pf ^ 
21° — 39° of stereoscopically triangulated loops with magnetic potential field 
models, and fiFFF ^ 15° — 21° for a force- free field model, which is partly 
caused by stereoscopic uncertainties fisE ^ 9°. We predict optimum conditions 
for solar stereoscopy during the time intervals of 2012-2014, 2016-2017, and 
2021-2023. 

Keywords: Sun: Corona — Stereoscopy 
1. Introduction 

Ferdinand Magellan's expedition was the first that completed the circumnavi- 
gation of our globe during 1519-1522, after discovering the Strait of Magellan 
between the Atlantic and Pacific ocean in search for a westward route to the 
"Spice Islands" (Indonesia), and thus gave us a first 360° view of our planet 
Earth. Five centuries later, NASA has sent two spacecraft of the STEREO 
mission on circumsolar orbits, which reached in 2011 vantage points on opposite 
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sides of the Sun that give us a first 360° view of our central star. Both discovery 
missions are of similar importance for geographic and heliographic charting, and 
the scientific results of both missions rely on geometric triangulation. 

The twin STEREO/A(head) and B(ehind) spacecraft (Kaiser et al. 2008), 
launched on 2006 October 26, started to separate at end of January 2007 by 
a lunar swingby and became injected into a heliocentric orbit, one propagating 
"ahead" and the other "behind" the Earth, increasing the spacecraft separation 
angle (measured from Sun center) progressively by about 45° per year. The two 
spacecraft reached the largest separation angle of 180° on 2011 February 6. A 
STEREO SECCHI CORl-A/B intercalibration was executed at 180° separation 
(Thompson et al. 2011). Thus, we are now in the possession of imaging data from 
the two STEREO/EUVI instruments (Howard et al. 2008; Wiilser et al. 2004) 
that cover the whole range from smallest to largest stereoscopic angles and 
can evaluate the entire angular range over which stereoscopic triangulation is 
feasible. It was anticipated that small angles in the order of ~ 10° should be 
most favorable, similar to the stereoscopic depth perception by eye, while large 
stereoscopic angles that are provided in the later phase of the mission would be 
more suitable for tomographic 3D reconstruction. 

The first stereoscopic triangulations using the STEREO spacecraft have been 
performed for coronal loops in active regions, observed on 2007 May 9 with 
a separation angle of asep = 7.3° (Aschwanden et al. 2008) and observed on 
2007 June 8 with asep = 12° (Feng et al. 2007). Further stereoscopic triangu- 
lations have been applied to oscillating loops observed on 2007 June 26 with a 
stereoscopic angle of asep = 15° (Aschwanden 2009), to polar plumes observed 
on 2007 Apr 7 with asep = 3.6° (Feng et al. 2009), to an erupting filament 
observed on 2007 May 19 with asep = 8.5° (Liewer et al. 2009), to an erupting 
prominence observed on 2007 May 9 with asep = 7.3° (Bemporad 2009), and 
to a rotating, erupting, quiescent polar crown prominence observed on 2007 
June 5-6 with asep = 11.4° (Thompson 2011). Thus, all published stereoscopic 
triangulations have been performed within a typical (small) stereoscopic angular 
range of asep « 3° — 15°, as it was available during the initial first months of 
the STEREO mission. The largest stereoscopic angle used for triangualtion of 
coronal loops was used for active region 10978, observed on 2007 December 
11, with a spacecraft separation of asep = 42.7° (Aschwanden and Sandman 
2010; Sandman and Aschwanden 2011), which produced results with similar 
accuracy as those obtained from smaller stereoscopic angles. So there exists 
also an intermediate rangle of aspect angles that can be used for stereoscopic 
triangulation. 

However, nothing is known whether stereoscopy is also feasible at large an- 
gles, say in the range of asep ^ 50° — 180°, and how the accuracy of 3D 
reconstruction depends on the aspect angle, in which range the stereoscopic 
correspondence problem is intractable, and whether stereoscopy at a maximum 
angle near asep ~ 180° is equally feasible as for asep ^ 0° for optically thin 
structures (as it is the case in soft X-ray and EUV wavelengths), due to the 
180° symmetry of line-of-sight intersections. In this study we are going to explore 
stereoscopic triangulation of coronal loops in the entire range of asep ^ 6° — 170° 
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Figure 1. Schematic figure of the spacecraft orbits of STEREO/A and B relative to Earth, 
with the spacecraft separation angles o-sep = ola—olb indicated approximately at the beginning 
of the years, ranging from 5° in April 2007 to ~ 180° in February 2011. 



and quantify the accuracy and quality of the results as a function of the aspect 
angle. 

Observations and data analysis are reported in Section 2, while a discussion 
of the results is given in Section 3, with conclusions in Section 4. 



2. OBSERVATIONS AND DATA ANALYSIS 

2.1. Observations 

We select STEREO observations at spacecraft separation angles with increments 
of ^ 45° over the range of asep ~ 0° to ~ 180°, which corresponds to time 
intervals of about a year during the past mission lifetime 2007-2011. A geometric 
sketch of the spacecraft positions STEREO/ A+B relative to the Earth-Sun axis 
is shown in Fig. 1. Additional constraints in the selection are: (i) The presence 
of a relatively large prominent active region; (ii) a position in the field-of-view of 
both spacecraft (since the mutual coverage overlap drops progressively from 180° 
initially to 0° during the first 4 years of the mission); (iii) a time near the central 
meridian passage of an active region viewed from Earth (to minimize confusion 
by foreshortening); and (iii) the availability of both STEREO/EUVI/A+B and 
calibrated SOHO/MDI data. The selection of 5 datasets is listed in Table 1, 
which includes the following active regions: (1) NOAA 10953 observed on 2007 
April 30 (also described in DeRosa et al. 2009; Sandman et al. 2009, Aschwanden 
and Sandman 2010; Sandman and Aschwanden 2011, Aschwanden et al. 2012), 
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Table 1. Data selection of 5 active regions observed with STEREO/EUVI and SOHO/MDI. 



Active 


Observing 


Observing 


Wave- 


Number of 


Spacecraft 


Region 


date 


times 


length 


stacked images 


separation 






(UT) 


(A) 


Na, Nb 


angle (deg) 


10953 (E23S10) 


2007-Apr-30 


23:00-23:20 


171 


8, 8 


6.1° 


10978 (E14S01) 


2007-Dec-ll 


16:30-16:50 


171 


8, 8 


42.7° 


11010 (E05N18) 


2009-Jan-12 


00:30-01:00 


171 


6, 6 


89.3° 


11032 (W03N16) 


2009-NOV-21 


00:30-01:00 


195 


12, 12 


126.9° 


11127 (W08N25) 


2010-NOV-23 


00:30-01:00 


195 


12, 6 


169.4° 



(2) NOAA region 10978 observed on 2007 December 11 (also described in As- 
chwanden and Sandman 2010, Aschwanden et al. 2012, and subject to an ongoing 
study by Alex Engell and Aad Van Ballegooijen, private communication), (3) 
NOAA 11010 observed on 2009 Jan 12, (4) NOAA 11032 observed on 2009 
Nov 21, and (5) NOAA 11127 observed on 2010 Nov 23. This selection covers 
spacecraft separation angles of agep ^ 6°, 43°, 89°, 127°, and 170°. 

For each of the 5 datasets we stacked the images during a time interval of 
20-30 minutes in order to increase the signal-to- noise ratio of the EUVI images. 
During the first 3 years (2006-2008) the nominal cadence of 171 A images was 
150 s, which yields 8 stacked images per 20 minute interval. Later in the mission, 
the highest cadence was chosen for the 195 A wavelength, but dropped from 150 
s to 300 s due to the reduced telemetry rate at larger spacecraft distances, which 
yields 6-12 stacked images per 30 minute interval. In one case (2010 Nov 23) 
the cadence in EUVI/ A and B are not equal, either due to data loss or different 
telemetry priorities (see time intervals and number of stacked images in Table 
1). The solar rotation during the time interval of stacked image sequences was 
removed to first order by shifting the images by an amount corresponding to the 
rotation rate at the extracted subimage centers. 

2.2. Stereoscopic Triangulation 

The geometric principles of stereoscopic triangulation are described in Aschwan- 
den et al. (2008; Sections 3.1 and 3.2 therein) for the general case of different 
spacecraft distances dA and ds from the Sun. The given formulas work correctly 
up to spacecraft separation angles of agep ~ 90°, but there is a sign ambiguity 
for larger separation angles. However, using the publicly available SSW/IDL 
software in the framework of the Wold Coordinate System (WCS) (Thompson 
2006), we can transform a pair of STEREO/ A-hB images into epipolar coordi- 
nates (Inhester 2006), by coaligning to the same Sun center position, derotating 
the spacecraft roll angles, and rescaling to the same solar distance d^), where 
stereoscopic triangulation is most straightforward. If (xa^Va) are the cartesian 
coordinates of a loop position in a STEREO/ A image (in units of solar radii 
measured from Sun center), the relationship between the heliocentric longi- 
tudes (Ia^Ib) and latitudes {bA^bs) in the images A and 5, and the cartesian 
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coordinates (xb^Vb) in the image B are (in the far-field approximation), 

Ba = arcsin (?/yi/r) 

Ia = arcsin [^^/(r cos6a)] 

Ib = + O^sep /-i \ 

bB =bA ' ^ ^ 

Xb = r sin(/B) cos(6b) 
Vb = VA 

where r = l-\-his the stereoscopically triangulated distance from Sun center (in 
units of solar radii). The conversion of image A pixel coordinates (m, Ja) into 
dimensionless Sun center coordinates (xa^Va) is, with Tpix the solar radius in 
units of EUVI pixels sizes, 



and the back- transformation into image B pixel coordinates is, 



where (io^jo) are the pixel coordinates of the Sun center in the co-registered 
epipolar STEREO image pair. There is no sign ambiguity in the coordinate 
transformation, as long as the triangulated feature is in front of the visible 
hemisphere (and plane-of-sky through Sun center) for each spacecraft, or more 
specifically, if the triangulated positions have positive values za > and zb > 
in each spacecraft coordinate system. 



2.3. Magnetic Field Models 



XA = {iA - io)/rpix 

VA = Ua -jo)/ Tpix ' 



iB = ^0 + Xb Tpix 
JB = jo + Vb Tpix ' 



(2) 
(3) 



Since the two STEREO/ A+B spacecraft have an almost symmetric separation 
angle in east and west direction with respect to Earth, their overlapping field- 
of-view is always centered closely to the Sun's central meridian for spacecraft 
separation angles of asep < 180°, and thus we have always also a magnetogram 
from an Earth-bound satellite available, such as from SOHO/MDI, during the 
considered time period of 2006-2011. We make use of the Michelson Doppler 
Imager (MDI) daily magnetic field synoptic full-disk data, which are taken every 
96 minutes, and thus are near-simultaneous with the STEREO/EUVI images 
within ^ 1 hour. The maximum magnetic field strengths of the 5 analyzed 
active regions are listed in Table 2, reaching up B ^ 3100 G. 

We calculate a magnetic potential-field model of each analyzed active region 
by the method of buried unipolar magnetic charges, which is described to first 
approximation in Aschwanden and Sandman (2010), and with a higher accuracy 
including the curvature of the solar surface in Aschwanden et al. (2012). Es- 
sentially, a line-of-sight magnetogram Bz{x,y) is decomposed into a number 
(typically ric = 200) components of buried magnetic charges, each one pa- 
rameterized with 4 parameters {Bj^xj^yi^ Zj) that characterize the surface field 
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strength Bj and the 3D position of the magnetic charge below the photospheric 
surface. The coronal potential field B(x) is calculated from the superposition of 
all Njn magnetic charges, 

B(x) = ^B,(x) = ^S,r^ (4) 

with r = \{x — Xj)^ {y — yj)^{z — zj)] being the distance of a coronal position 
[x^y^z] from the magnetic charge j, vj = (xj^yj^Zj) is the subphotospheric 

position of the buried charge, and dj = ^1 — x'j — y'j — Zj is the depth of the 

magnetic charge. 

We calculate magnetic field lines by starting at photospheric footpoints and 
extrapolating along the local magnetic field vector B(x) = {B^^ By^Bz) in steps 
of As = 0.004 solar radii. We initiate the field line footpoints in a regular grid 
(of say 50 x 50) footpoints positions, but plot only those field lines that have a 
photospheric field strength above some threshold, say B > 100 G. 

For the first three active regions where we have sufficient constraints by 
stereoscopic loops {Nioops ~ 20), we calculate also a nonlinear force- free field 
(NLFFF) solution according to a new code based on an analytical approxima- 
tion of divergence-free and force-free fields that includes an azimuthal magnetic 
field component with vertical twist and is accurate to second-order (of the 
force- free parameter a). This new NLFFF code uses the constraints of a line- 
for-sight magnetogram to define the potential field (Eq. 4) and is suitable for 
fast forward- fitting to coronal field constraints, such as the stereoscopically tri- 
angulated 3D loop coordinates as calculated here. The analytical theory of this 
new NLFFF code is described in Aschwanden (2012), the numerical code and 
tests in Aschwanden and Malanushenko (2012), and first solar applications to 
four stereoscopically observed active regions in Aschwanden et al. (2012). 

2.4. Data Analysis and Results 

In Figs. 2-6 we present the results of stereoscopically triangulated loops for the 
5 analyzed active regions and compare them with magnetic potential field and 
force- free field models. In each of the Figures we use the following layout of pan- 
els: Partial STEREO/ A and B images are shown in the top panels, each one with 
a field-of-view that encompasses the active region of interest. The STEREO/ A 
and B images are shown for the brightness on a logarithmic color scale (top 
panels), as well as in a highpass-filtered version (by unsharp masking, i.e., by 
subtracting a 3 x 3 boxcar-smoothed image from the original image) to enhance 
the loop structures (panels in second row). The SOHO/MDI magnetogram is 
shown in the bottom right panel, which has a different field-of-view from an 
Earth-bound vantage point. In addition we show a side view of the active region 
by rotating the magnetogram view by 90° to the north (bottom left panel). In all 
panels, the projections of stereoscopically triangulated loops are shown in blue 
color, and the magnetic potential field lines in red color. 
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Table 2. Data analysis of selection 5 active regions. 



Active 


Number 


Misalignment 


Misalignment 


Stereoscopy 


Magnetic 


Region 


of 


angle 


angle 


error 


field strength 




loops 


/iPF (deg) 


liFFF (deg) 


IJ^SE (deg) 


B (Gauss) 


10953 


100 


27.7° 


19.8° 


9.4° 


[-3134,+1425] 


10978 


52 


20.8° 


15.1° 


8.9° 


[-2270,+2037] 


11007 


20 


36.2° 


16.5° 




[ -737,+1342] 


11035 


15 


29.8° 


17.4° 




[ -997, +774] 


11161 


5 


39.4° 


21.2° 




[-2033, +884] 



The first case is active region NOAA 10953 (Fig. 2), where we display the 
same 100 loop segments that have been triangulated in an earlier study (As- 
chwanden and Sandman 2010). The spacecraft separation angle is asep = 6.1° 
and the almost identical direction of the line-of-sights of both STEREO / A and B 
spacecraft makes it easy to identify the corresponding loops in A and B, and thus 
the triangulation is very reliable. Note that the height range where discernable 
loops can be traced in the highpass-filtered images is about hmax ^0.1 solar 
radii (or ^ 70 IVIm), which is commensurable with the hydrostatic density scale 
height expected for a temperature of T = 1.0 MK that corresponds to the peak 
sensitivity of the EUVI 171 A filter. This is particularly well seen in the side 
view shown in the bottom left panel in Fig. 2. A measurement of the mean 
misalignment angle averaged over 10 positions of the 100 reconstructed loops 
with the local magnetic potential field shows a value of /j^pf = 27.7° (Table 
2), similar to earlier work (Aschwanden and Sandman 2010; Sandman and As- 
chwanden 2011). However, forward- fitting of a nonlinear force- free field model 
reduces the misalignment to (Ifff = 19.8°, which implies that this active region 
is slightly nonpotential. The remaining misalignment is attributed to at least two 
reasons, partially to inadequate parameterization of the force-free field model, 
and partially to stereoscopic measurement errors asE due to misidentified loop 
correspondences and limited spatial resolution. An empirical estimate of the 
stereoscopic error was devised in Aschwanden and Sandman (2010), based on 
the statistical non-parallelity of closely-spaced triangulated loop 3D trajectories, 
which yielded for this case a value oi asE = 9.4°. In summary, we find that this 
active region is very suitable for stereoscopy, allows to discern a large number 
(100) of loops, minimizes the stereoscopic correspondence problem due to the 
small {asep = 6.1°) spacecraft separation angle, displays a moderate misalign- 
ment angle and stereoscopic measurement error {fisE = 9.4°). This well-defined 
case will serve as a reference for stereoscopy at larger angles. 

The second case is active region NOAA 10978 (Fig. 3), observed on 2007 
Dec 11 with a spacecraft separation angle of agep = 42.7°. Note that the views 
from EUVI/ A and B appear already to be significantly different with regard 
to the orientation of the triangulated loops, as seen from a distinctly different 
aspect angle. A set of 52 coronal loops were stereoscopically triangulated in 
this region (Aschwanden and Sandman 2010), a mean misalignment angle of 
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Figure 2. STEREO/EUVI spacecraft A (top right) and spacecraft B image (top left), of 
active region AR 10953, observed on 2007 Apr 30, 23:00 UT in the 171 A wavelength, with a 
spacecraft separation angle of asep = 6.1°. The images are highpass-filtered to enhance loop 
structures (middle left and right panels). A near-simultaneous SOHO/MDI magnetogram is 
shown (bottom right), overlaid with the stereoscopically triangulated loops (blue curves) and 
magnetic field lines computed with a nonlinear force- free model (red curves), viewed from the 
direction of Earth or SOHO/MDI (bottom right), and rotated by 90*-* to the north (bottom 
left). 
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Figure 3. Active region AR 10978, observed on 2007 Dec 11, 16:30 UT in the 171 A wave- 
length, with a spacecraft separation angle of a.sep = 42.7°. A SOHO/MDI magnetogram is 
shown (bottom right), overlaid with the stereoscopically triangulated loops (blue curves) and 
magnetic field lines computed with a nonlinear force- free model (red curves). 
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Figure 4. Active region AR 11010, observed on 2009 Jan 12, 00:30 UT in the 171 A wave- 
length, with a spacecraft separation angle of asep = 89.3°. A SOHO/MDI magnetogram is 
shown (bottom right), overlaid with the stereoscopically triangulated loops (blue curves) and 
magnetic field lines computed with a nonlinear force- free model (red curves). 
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Figure 5. Active region AR 11032, observed on 2009 Nov 21, 00:30 UT in the 195 A wave- 
length, with a spacecraft separation angle of asep = 126.9°. A SOHO/MDI magnetogram is 
shown, overlaid with the stereoscopically triangulated loops (blue curves) and magnetic field 
lines computed with a potential field model. Magnetic field lines have a footpoint threshold of 
B > 100 G. 
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Figure 6. Active region AR 11127, observed on 2010 Nov 23, 00:30 UT in the 195 A wave- 
length, with a spacecraft separation angle of asep = 169.4°. A SOHO/MDI magnetogram is 
shown, overlaid with the stereoscopically triangulated loops (blue curves) and magnetic field 
lines computed with a potential field model. Magnetic field lines have a footpoint threshold of 
B > 100 G. 
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lipF = 20.8° is found for a potential field model, and a reduced value of (Ifff = 
15.1° is found for the force- free model (Table 2), while a stereoscopic error of 
otsE = 8-9° is estimated (Aschwanden and Sandman 2010). Thus, the quality of 
stereoscopic triangulation (as well as the degree of non-potentiality) is similar to 
the first active region, although we performed stereoscopy with a 7 times larger 
spacecraft separation angle {asep =42.7°) than before {asep = 6.1°). Apparently, 
stereoscopy is still easy at such angles, partially helped by the fact that the active 
region is located near the central meridian (±20°) for both spacecraft, which 
provides an unobstructed view from top down, so that the peripheral loops of 
the active region do not overarch the core of the active region, where the bright 
reticulated moss pattern (Berger et al. 1999) makes it almost impossible to 
discern faint loops in the highpass-filtered images. The top-down view provides 
also an optimum aspect angle to disentangle closely-spaced loops, which is an 
important criterion in the stereosopic correspondence identification. 

The third case is active region NOAA 11010 (Fig. 4), observed on 2009 Jan 
12 with a near-orthogonal spacecraft separation angle of asep = 89.3°. Due 
to the quadrature of the spacecraft, only a sector of ±45° east and west of 
the central meridian (viewed from Earth) is jointly visible by both spacecraft. 
This particular active region is seen at a 45° angle by both STEREO/ A and 
STEREO/B. This symmetric view is the optimum condition to discern a large 
number of inclined loop segments and to identify the stereoscopic correspon- 
dence. We triangulate some 20 loop segments, which appear almost mirrored in 
the STEREO/ A and B image due to the east- west symmetry of the magnetic 
dipole. A mean misalignment angle oi jipF = 36.2° with the potential field model 
is found, and a reduced value oi iifff = 16.5° with the force- free field model. An 
estimate of the statistical (non-parallelity) stereoscopic error is not possible due 
to the small number of triangulated loops. Thus, we conclude that stereoscopy is 
still possible in quadrature. Mathematically, the orthogonal projections should 
yield the most accurate 3D coordinates of a curvi-linear structure, but in prac- 
tice, confusion of multiple structures with near-aligned projections can cause a 
disentangling problem in the stereoscopic correspondence identification at this 
intermediate angle. 

The fourth case is active region NOAA 11032 (Fig. 5), observed on 2009 Nov 
21 with a large spacecraft separation angle of agep = 126.9°. STEREO/ A sees 
the active region near the east limb from an almost side-on perspective, while 
STEREO/B sees a similar mirror image near the west limb, where confusion near 
the limb makes the stereoscopic correspondence identification more difficult. We 
trace some 15 loop segments, but do not succeed in pinning down a larger number 
of loops, partially because this active region is small and does not exhibit numer- 
ous bright loops, and partially because of increasing confusion problems near the 
limb. We searched for larger active regions over several months around this time, 
but were not successful due to a dearth of solar activity during this time. We find 
a misalignment angle of fipF = 29.8° for the potential field, and ijlfff = 17.4° 
for the force-free field model, which is still comparable with the previous active 
regions triangulated at smaller sterescopic angles. Thus, stereoscopy seems to be 
still feasible at such large stereoscopic angles. 
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The last case is active region NOAA 11127 (Fig. 6), observed on 2010 Nov 23 
with a very large spacecraft separation angle of agep = 169.4°, only two months 
before the two STEREO spacecraft pass the largest separation point. At this 
point, the common field-of-view that is overlapping from STEREO /A and B 
is only the central meridian zone seen from Earth (or the opposite meridian 
behind the Sun). STEREO/A observes active region NOAA 11127 at its east 
limb, while STEREO/B sees it at its west limb, so both spacecraft see only 
the vertical structure of the active region from a side view (see Fig. 6 top). 
This particular configuration is very unfavorable for stereoscopy. Although the 
vertical structure in altitude can be measured very accurately, the uncertainty 
in horizontal direction in longitude is very large and suffers moreover the sign 
ambiguity of positive or negative longitude difference with respect to the limb 
seen from Earth. Consequently, we have reliable information on the altitude and 
latitude of loops, while the longitude is essentially ill-defined. In order to reduce 
the large scatter in the measurement of ^-coordinates along a loop, introduced 
by the near-infinite amplification of parallax uncertainties tangentially at the 
limb, we restrict the general solution of geometric 3D triangulation to planar 
loops, by applying a linear regression fit of the z{y) coordinates. The example in 
Fig. 6 shows that we can trace some (5) loops in the plane of the sky and have 
no problem in identifying the stereoscopic counterparts in both STEREO/ A 
and B images, but the stereoscopic triangulation is ill-defined at this singularity 
of the sign change in the parallax effect. The misalignment between the three 
loop directions and the potential field is jj^pf = 39.4°, and for the force- free 
field model is /hfff = 21.2°, which indicates that the orientation of the loop 
planes is less reliably determined. Stereoscopic triangulation brakes down at 
this singularity of separation angles at agep ^ 180°, although the stereoscopic 
correspondence problem is very much reduced for the "mirror images", similar 
to the near-identical images at small separation angles asep ~ 0°. 



3. DISCUSSION 

We are discussing now the pro's and con's of stereoscopy at small and large aspect 
angles, which includes quantitative estimates of the formal error of stereoscopic 
triangulation (Section 3.1), the stereoscopic correspondence and confusion prob- 
lem (Section 3.2), and the statistical probability of stereoscopable active regions 
during the full duration of the STEREO mission (Section 3.3), all as a function 
of the stereoscopic aspect angle (or spacecraft separation angle asep in the case 
of the STEREO mission). 

3.1. Stereoscopic Triangulation Error 

Stereoscopic triangulation involves a parallax angle around the normal of the 
epipolar plane. For the STEREO mission, the epipolar plane intersects the Sun 
center and the two spacecraft A and B positions, which are separated mostly 
in east-west direction. No parallax effect occurs when the loop axis coincides 
with the epipolar plane, i.e., when the loop axis points in east- west direction. 
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Thus, the accuracy of stereoscopic triangulation depends most sensitively on this 
orientation angle i^^ which we define as the angle between the loop direction and 
the normal of the epipolar plane (i.e., approximately the y-axis of a solar image 
in north-south direction). If the position of a loop centroid can be determined 
with an accuracy of a half pixel size the dependence of the stereoscopic 

error on the orientation angle is then (Aschwanden et al. 2008), 

Thus, for a highly inclined loop that has an eastern and western foot point at 
the same latitude, the stereoscopic error is minimal near the loop footpoints 
(pointing in north-south direction) and at maximum near the loop apex (pointing 
in east- west direction). 

In addition, the accuracy of stereoscopic triangulation depends also on the 
aspect angle (or spacecraft separation angle a^ep), which can be quantified with 
an error trapezoid as shown in Fig. 7. If the half separation angle agepf^ is defined 
symmetrically to the Earth-Sun axis (z-axis), the uncertainty in z-direction is 
az — {/S.pix/2) / ^m.{asep/'2) and is in the x-direction (in the epipolar plane) 
(Jx = (Ap^a;/2)/ sin([7r — asep]/2). Including also a half pixel-size error ay = 
{Apix/2) in the y-direction, we have then a combined error for the 3D position 
of a triangulated point as, 

cr{asep) = \ cfI 

z y sm-[^asep/'^) sm" [^tt - a^epj/^^J 

(6) 

This positional error is symmetric for small asep and large stereoscopic angles 
(tt — OLsep)^ and has a minimum at an orthogonal angle of asep = 7r/2. The error 
is largest in z-direction for small spacecraft separation angles, while it is largest 
in x-direction for separation angles near agep ~ 180°. 

To compare the relative importance of the two discussed sources of errors we 
can evaluate the parameters that increase the individual errors by a factor of 
two. This is obtained when the orientation angle of a loop segment (with respect 
to the east- west direction) increases from ^ = {)° to 60°, or if the spacecraft 
separation angle changes from the optimum angle asep = 90° to asep = 25° 
(or asep = 155°, respectively). If stereoscopy at small angles of asep = 5° is 
attempted, the positional error is about ten- fold (corresponding to ^ SApix), 
compared with the optimum angle at asep = 90° (corresponding to ^ O.SApix). 
For STEREO/EUVI with a pixel size of Apix = 1.6'' ^1.2 Mm, this amounts 
to an accuracy range of ^ 1 — 10 Mm. 

3.2. Stereoscopic Correspondence and Confusion Error 

The previous considerations are valid for isolated loops that can be unambigu- 
ously disentangled in an active region, in both the STEREO/ A and B images. 
However, this is rarely the case. In crowded parts of active regions, the corre- 
spondence of a particular loop in image A with the identical loop in image B 




SOLA: ms.tex; 13 July 2012; 0:15; p. 15 



M.J. Aschwanden et al. 




A 



Az 



Figure 7. The error trapezoid of stereoscopic triangulation is shown (grey area), given by 
the two hne-of-sights of the two observer directions A and B, separated by an angle asep- The 
uncertainties Ax in x-direction and spacecraft in z-direction depend on the pixel width Apix 
and half aspect angle asep/'^- 

can often not properly be identified. We call this confusion problem also the 
stereoscopic correspondence problem^ which appears in every stereoscopic tie- 
point triangulation method. In order to quantify this source of error, we have to 
consider the area density of loops and their relative orientation. A top-down view 
of an active region, e.g., as seen for small stereoscopic angles by both spacecraft 
for an active region near disk center (e.g.. Fig. 2), generally allows a better 
separation of individual loops, because only the lowest density scale height is 
detected (due to hydrostatic gravitational stratification), and neighbored loop 
segments do not obstruct each other due to the foreshortening projection effect 
near the footpoints. In contrast, every active region seen near the limb, shows 
many loops at different longitudes, but at similar latitudes, cospatially on top 
of each other, which represents the most severe confusion problem. Thus, we 
can essentially quantify the degree of confusion by the loop number density per 
pixel, which approximately scales with the inverse cosine-function of the center- 
to-limb angle due to foreshortening. In other words, we can define a quality factor 
C{p) for identifying the stereoscopic correspondence C of properly disentangled 
coronal loops, which drops from C(p = 0) = 1 at disk center to C(/? = 7r/2) = 
at the limb, where p is the center-to-limb angle measured from Sun center. 



The orbits of the STEREO mission reduce the overlapping area on the solar 
surface that can be jointly viewed by both spacecraft A and B linearly with 
increasing spacecraft separation angle, so that the center-to-limb distance p of 



C{p) = \cos{p)\ 



0<\p\< 7r/2 . 



(7) 
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Figure 8. Quality factor Q of stereoscopic triangulation as a function of the spacecraft sepa- 
ration angle asep^ which is a function of the accuracy A of triangulated stereoscopic positions 
and the stereoscopic correspondence quality factor C. The best quality (within a factor of 2) 
occurs in the range of a.sep = 22° — 125°. 



an active region located on the central meridian (viewed from Earth) is related 
to the spacecraft separation angle by 

P=^, (8) 



increasing linearly with the separation angle from p = at the beginning of the 
mission to p = 7r/2 at maximum spacecraft separation angle. The location of an 
active region at the central meridian provides the best view for both spacecraft, 
because an asymmetric location would move the active region closer to the limb 
for one of the spacecraft, and thus would increase the degree of confusion, as we 
verified by triangulating a number of asymmetric cases. Thus, we can express 
the quality factor of stereoscopic correspondence C (Eq. 7) by the spacecraft 
separation angle (Eq. 8) and obtain the relationship, 

C{asep) = cos(asep/2) < asep < TT . (9) 

Defining a quality factor Q for stereoscopic triangulation by combining the 
stereoscopic correspondence quality C (Eq. 9) with the accuracy A of stereo- 
scopic positions (Eq. 6), which we may define by the normalized inverse error 
(i.e., A = cF{asep,min)/cr{asep)), we obtain 

=C.A= ^^-o^M2) ^^^^ 

y^l + 1/ sin^ (a,ep/2) + 1/ sin^ [(tt - a,ep)/2] 
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Figure 9. The statistically expected variation of the number of active regions N(t) during 
the solar cycle (thin solid curve), the spacecraft separation angle asep(t) (dotted curve), the 
stereoscopic accuracy A(t) (dashed curve), and the expected number of stereoscopable active 
regions N^Rit) (curve with grey area) as a function of time during a full 16-year mission cycle 
of the STEREO mission. 



We plot the functional dependence of this stereoscopic quality factor Q{asep) 
together with their underlying factors A{asep) and C{asep) in Fig. 8 and obtain 
now an asymmetric function of time (or spacecraft separation angle) that favors 
smaller stereoscopic angles. The stereoscopic quality factor is most favorable 
(within a factor of 2) in the range of agep = 22° — 125°, which corresponds to 
the mission phase between August 2007 and November 2009. The same optimum 
range will repeat again at the backside of the Sun 5 years later between August 
2012 and November 2014. 

From our analysis of 5 active regions spread over the entire spacecraft separa- 
tion angle range we find acceptable results regarding triangulation accuracy in 
the range of separation angles of agep = 6° and 127° (based on acceptable 
misalignment angles of apF ^ 35°, which coincides with the predicted opti- 
mum range of agep ^ 22° — 125°, while stereoscopy definitely brakes down at 
asep ~ 170°, as predicted by theory (Fig. 8 and Eq. 10). 

3.3. The Statistical Probability of Stereoscopy 

There are different factors that affect the quality or feasibility of solar stere- 
oscopy, such as (i) the availability of large active regions (which varies sta- 
tistically as a function of solar cycle), (ii) the simultaneous viewing by both 
spacecraft STEREO/ A and B (which depends on the spacecraft separation 
angle), (iii) the geometric foreshortening that affects the stereoscopic correspon- 
dence problem (which depends on the center-to-limb distance for each spacecraft 
view), and (iv) the time of the central meridian passage of the active region 
for a viewer from Earth (which determines the symmetry of views for both 
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spacecraft, where minimum confusion occurs in the stereoscopic correspondence 
identification). Ah but the first factor depend on the spacecraft separation angle 
<^sep{t)^ which is a specific function of time for the STEREO mission (with a 
complete cycle of 16 years). In order to assess the science return of the STEREO 
mission or future missions with stereoscopic capabilities, it is instructive to quan- 
tify the statistical probability of acceptable stereoscopic results as a function of 
spacecraft separation angle or time. 

We already quantified the quality of stereoscopy Q{asep) as a function of the 
stereoscopic angle in Eq. (10). Let us define the number probability N{t) of 
existing active regions at a given time t with a squared sinusoidal modulation 
during the solar cycle. 



where A^o ^ 10 is the maximum number of active regions existing on the total 
solar surface during the maximum of the solar cycle, to is the time of the solar 
minimum (e.g., tmin ^ 2009), and Tcyde ^ 10.3 yrs the current average solar 
cycle length. 

The second effect is the overlapping area on the solar surface that is simul- 
taneously seen by both spacecraft STEREO/ A and B, which decreases linearly 
with the spacecraft separation angle from 50% at asep{t = ti) (with ti = 2007.1 
at the start of spacecraft separation) to 0% at asep 

{t = t2) (with t2 = 2011.1 
at maximum separation), and then increases linearly again for the next quarter 
phase of a mission cycle. If we fold the variation N(t) of the solar cycle (Eq. 11) 
with the triangular stereoscopic overlap area variation A{t) together, we obtain 
a statistical probability for the number of stereoscopically triangulable active 
regions. However, the number of accurate stereoscopical triangulations scales 
with the quality factor Q{asep) (Eq. 10), where the spacecraft separation angle 
asep(t) is a piece- wise linear (triangular) function of time t according to the 
spacecraft orbit. Essentially we are assuming that the probability of successfull 
stereoscopic triangulations at a given time scales with the quality factor or feasi- 
bility of accurate stereoscopy at this time. So, we obtain a combined probability 
of stereoscopically triangulable active regions of. 



In Fig. 9 we show this combined statistical probability of feasible stereoscopy 
in terms of the expected number of active regions for a full mission cycle of 16 
years, from 2006 to 2022. It shows that the best periods for solar stereoscopy 
are during 2012-2014, 2016-2017, and 2021-2023. 

4. CONCLUSIONS 

After the STEREO mission reached for the first time a full 360° view of the Sun 
this year (2007 Feb 6), the two STEREO A and B spacecraft covered also for the 
first time the complete range of stereoscopic viewing angles from agep ~ 0° to 




(11) 



NAR{t) = N{t) X Q[asep{t)] . 



(12) 
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asep ~ 180°. We explored the feasibility of stereoscopic triangulation for coronal 
loops in the entire angular range by selecting 5 datasets with viewing angles at 
asep ^ 6°, 43°, 89°, 127° and 169°. Because previous efforts for solar stereoscopy 
covered only a range of small stereoscopic angles {asep ~ 45°), we had to gener- 
alize the stereoscopic triangulation code for large angles up to asep < 180°. We 
find that stereoscopy of coronal loops is feasible with good accuracy for cases 
in the range asep ~ 6° — 127°, a range that is also theoretically predicted by 
taking into account the triangulation errors due to finite spatial resolution and 
confusion in the stereoscopic correspondence identification in image pairs, which 
is hampered by projection effects and foreshortening for viewing angles near the 
limb. Accurate stereoscopy (within a factor of 2 of the best possible accuracy) is 
predicted for a spacecraft separation angle range of asep ^ 22° — 125°. Based on 
this model we predict that the best periods for stereoscopic 3D reconstruction 
during a full 16-year STREREO mission cycle occur during 2012-2014, 2016- 
2017, and 2021-2023, taking the variation in the number of active regions during 
the solar cycle into account also. 

Why is the accuracy of stereoscopic 3D reconstruction so important? Solar 
stereoscopy has the potential to quantify the coronal magnetic field indepen- 
dently of conventional 2D magnetogram and 3D vector magnetograph extrap- 
olation methods, and thus serves as an important arbiter in testing theoretical 
models of magnetic potential fields, linear force- free field models (LFFF), and 
nonlinear force- free field models (NLFFF). A benchmark test of a dozen of 
NLFFF codes has been compared with stereoscopic 3D reconstruction of coronal 
loops and a mismatch in the 3D misalignment angle of /i ^ 24° — 44° has 
been identified (DeRosa et al. 2009), which is attributed partially to the non- 
force-freeness of the photospheric magnetic field, and partially to insufficient 
constraints of the boundary conditions of the extrapolation codes. Empirical 
estimates of the error of stereoscopic triangulation based on the non-parallelity 
of loops in close proximity has yielded uncertainties of fisE ^ 7° — 12°. Thus the 
residual difference in the misalignment is attributed to either the non-potentiality 
of the magnetic field (in the case of potential field models), or to the non- 
force- freeness of the photospheric field (for NLFFF models). We calculated also 
magnetic potential fields here for all stereoscopically triangulated active regions 
and found mean misalignment angles of /ipi? ~ 21° — 39°, which improved to 
fiFFF 15° — 21° for a nonlinear force- free model, which testifies the reliability 
of stereoscopic reconstruction for the first time over a large angular range. The 
only case where stereoscopy clearly fails is found for an extremely large sepa- 
ration angle of {asep ^ 170°), which is also reflected in the largest deviation of 
misalignment angles found {jiNP ^ 39°, (Ifff ^ 21°). Based on these positive 
results of stereoscopic accuracy over an extended angular range from small to 
large spacecraft separation angles we anticipate that 3D reconstruction of coronal 
loops by stereoscopic triangulation will continue to play an important role in 
testing theoretical magnetic field models for the future phases of the STEREO 
mission, especially since stereoscopy of a single image pair does not require a 
high cadence and telemetry rate at large distances behind the Sun. 
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